Thermal radiation and amplified spontaneous emission from a random medium 
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We compute the statistics of thermal emission from systems in which the radiation is scattered 
chaotically, by relating the photocount distribution to the scattering matrix — whose statistical 
properties are known from random-matrix theory. We find that the super-Poissonian noise is that 
of a black body with a reduced number of degrees of freedom. The general theory is applied to a 
disordered slab and to a chaotic cavity, and is extended to include amplifying as well as absorbing 
systems. We predict an excess noise of amplified spontaneous emission in a random laser below the 
laser threshold. 
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The emission of photons by matter in thermal equilib- 
rium is not a series of independent events. The text- 
book example is black-body radiation Consider 
a system in thermal equilibrium (temperature T) that 
fully absorbs any incident radiation in N{lo) propagat- 
ing modes within a frequency interval Slo around lj. A 
photodetector counts the emission of n photons in this 
frequency interval during a long time t 3> \/5lo. The 
probability distribution P{n) is given by the negative- 
binomial distribution with v — NtSLo/^n degrees of free- 
dom, 



P(ri) cx 



cxp{—nhuj / k^T) . 



(1) 



The binomial coefficient counts the number of partitions 
of n bosons among v states. The mean photocount 
n = vf is proportional to the Bose-Einstein function 



(2) 



In the limit fi/v 0, Eq. ([^) approaches the Poisson 
distribution P{n) cx n"" /n\ of independent photocounts. 
The Poisson distribution has variance Var n = n equal to 
its mean. The negative-binomial distribution describes 
photocounts that occur in "bunches" , leading to an in- 
crease of the variance by a factor 1 + n/iy. These basic 
facts are known since the beginning of this century jsj . 

Thermal radiation is also referred to as "chaotic ra- 
diation" In recent years the word "chaotic" has 
entered optics in a different context, to describe systems 
that scatter radiation in an irregular, random way [Q. 
Such systems, typically, have weak absorption, so they 
are far from being black bodies. Two recent papers have 
studied deviations from black-body radiation in the case 
of one-dimensional scattering H,^, but chaotic systems 
are intrinsically not one-dimensional. What, then, is the 
statistics of the chaotic radiation resulting from chaotic 
scattering? That is the problem addressed in this paper. 

This problem is significant for more than one reason. 
First of all, thermal emission is a fundamental property 
of a system. Deviations from the black-body limit con- 
tain information on chaotic scattering that can not be 



obtained from classical scattering experiments. Most 
studies of the optical properties of random media have 
been restricted to classical optics [Q. The similarity be- 
tween the classical wave equation and the Schrodinger 
equation has permitted the transfer to classical optics 
of powerful theoretical techniques from condensed mat- 
ter physics Our solution of the thermal-radiation 
problem demonstrates how one of these techniques, the 
method of random-matrix theory can be applied to 
quantum optics. That is the second reason for the sig- 
nificance of this problem. The third reason is the recent 
interest in amplifying random media, motivated by pos- 
sible applications as a "random laser" ||l^,|ll|. A linear 
amplifier can be thought of as being in thermal equilib- 
rium at a negative temperature p^ , so that our theory 
of thermal radiation can also deal with amplified sponta- 
neous emission. 




Fig. 1^ — Schematic diagram of a random medium (dot- 
ted) connected to a photodetector (shaded) via an A^- 
mode waveguide. 

We start with the formulation and solution of the prob- 
lem in general form, and then turn to specific appli- 
cations. We consider a random medium coupled to a 
photodetector via a waveguide (in vacuum) with N{uj) 
propagating modes (counting polarizations) at frequency 
uj. (See Fig. |l|.) We assume that any Brownian motion 
of the scattering centra in the random medium can be 
disregarded on the time scale of the measurements. The 
scattering rate is denoted by 1/ts, and the absorption or 
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amplification rate by l/ra- To quantize the electromag- 
netic field we use the method of input-output relations 
developed by Gruner and Welsch |^ and by Loudon and 
coworkers |^,^,|3|. The incoming and outgoing modes 
in the waveguide are represented by two A^-component 
vectors of annihilation operators a"^{uj), a°"'(a;). They 
satisfy the commutation relations 

[a„(w), aj„(w')] ^ Snm5{uj - uj'), [a„(w), am(w')] =0, 



for 



a or a = a 



(3) 

The input-output relations take 



the form |g,|Tl,|ll 

^out = s -a"' + U -b+V -c^ 



(4) 



with S{llj) the N x N scattering matrix. The boson op- 
erators b and c satisfy Eq. (pi) provided 



(5) 



(1 denoting the iVxA^ unit matrix). The matrix 1 — 5'-S'''^ 
is positive definite in an absorbing medium, so we can put 

= 0. Conversely, in an amplifying medium i — S ■ is 
negative definite, so we can put U — 0. This determines 
U, y up to a unitary transformation. All our final expres- 
sions depend only on the combination U ■ U''' — V ■ V\ so 
that any freedom in the choice of U, V is irrelevant once 
the scattering matrix is fixed. 

Eq. (H) can be understood as a fluctuation-dissipation 
relation: The left-hand side accounts for quantum fluc- 
tuations in the electromagnetic field due to spontaneous 
emission or absorption of photons, the right-hand side 
accounts for dissipation due to absorption (or stimulated 
emission in the case of an amplifying medium). Eq. (|^) 
also represents a link between classical optics (the scat- 
tering matrix S) and quantum optics (the quantum fluc- 
tuation matrices U, V) . 

In an absorbing medium, the operator b accounts for 
thermal emission with expectation value 



{bi{u;)b^{io'))=6nraS{uj-Lu')f{u;,T). 



(6) 



The inverted oscillator c accounts for spontaneous emis- 
sion in an amplifying medium. We consider the regime of 
linear amplification, below the laser threshold. Formally, 
this regime can be described by a thermal distribution at 
negative temperature — T, 



(C„(w)cj„(w')) -5n,nS{uJ - Uj')f{uJ, -T), 



(7) 



the zero-temperature limit corresponding to a complete 
population inversion [T^ . Higher order expectation val- 
ues are obtained by pairwise averaging, as one would do 
for Gaussian variables, after having brought the opera- 
tors into normal order. 

The incoming radiation is in the vacuum state, while 
the outgoing radiation is collected by a photodetector 
[Q. The probability that n photons are counted in a 
time t is given by |l5| , [T^ 



P{n) 



re 



(8) 



(The colons denote normal ordering.) It is convenient to 
work with the generating function F(^) = Kp^P/pl of 
the factorial cumulants k„ 11711, 



(9) 



To evaluate F{^) we substitute Eq. (|[) into Eq. (||) and 
perform the Gaussian averages. 

A simple expression results in the long-time regime, 

/•OC 1 

F{0=-t -^^ln||l-(l-5-5t)C/||, (10) 
Jq ^tt 

where || • ■ ■ || indicates the determinant. Eq. (^) is valid 
when LUct 3> 1, with lUc the frequency interval within 
which S ■ iS^ does not vary appreciably. We have also 
found a simple expression in the short-time regime, 

/•OO 7 

F{0 = -\n\\l-tJ^ ^{1-S-S^)^fl (11) 

valid when flf-t <ti 1, with fi^ the frequency range over 
which S ■ S*^ differs appreciably from the unit matrix. 
(The reciprocal of Oc is the coherence time of the thermal 
emissions.) The two equations (|l^) and (^ij) are the key 
results of this paper. They reduce the quantum optical 
problem of the photon statistics to a computation of the 
scattering matrix of the classical wave equation. That is 
a major simplification, because the statistical properties 
of the scattering matrix of a random medium are known 
from random-matrix theory ^8 19 1. 

The long-time limit (^o|) is particularly simple, as it 
depends only on the set of eigenvalues ai, a2, ■ ■ ■ o'n of 
S ■ . We call the cr„'s "scattering strengths". An 
additional simplification of the long-time regime is that 
one can do a frequency- resolved measurement, counting 
only photons within a narrow frequency interval Sio (with 
LUc ^ ^ 1 /t) ■ The factorial cumulants are then given 
by 



JV 



^p = (p-l)!^.r7V-l^(l-a„)^ 



(12) 



where v = Nt5ui/2TT was defined in the introduction. 
For comparison with black-body radiation we parame- 
terize the variance in terms of the effective number 
of degrees of freedom 0] , 



VarTi = n{l + n/vcs), 



(13) 
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with VcS = V for a black body. Eq. ( |T^ implies 



< 1. 



(14) 



We conclude that the super-Poissonian noise of a ran- 
dom medium corresponds to a black body with a reduced 
number of degrees of freedom. Note that the reduction 
occurs only for > 1. 

We now turn to applications of our general formulas to 
specific random media. We concentrate on the long-time, 
frequency- resolved regime with ^ 1, leaving the short- 
time and single-mode regimes, and the case of broad- 
band detection, for future publication |^^. An ensemble 
of random media has a certain scattering-strength den- 
sity p(cr). For A^ ^ 1 sample-to-sample fluctuations are 
small, so the ensemble average is representative for a sin- 
gle system. We nmy therefore replace ^„ by J da p{a) 
in Eqs. (p|) and (|i? 




Fig. 1^ — Effective number of degrees of freedom as a 
function of normalized absorption or amplification rate. 
The dashed curve is for the disordered slab, the solid 
curves are for the chaotic cavity. The amplifying slab 
would be above the laser threshold for any 7, so we only 
plot the case of absorption. For the cavity both the cases 
of absorption and amplification are shown. The black- 
body limit for absorbing systems and the laser threshold 
for amplifying systems are indicated by arrows. 

As a first example we compute the thermal radiation 
from a disordered absorbing slab. The slab is sufhciently 
thick that there is no transmission through it, represent- 
ing a semi-infinite random medium. We define the nor- 
mahzed absorption rate ||2|] 7 — ^Tg/ra. The scattering- 
strength density p(a) in the regime jN'^ ^ 1 is known 
1^,^. It is non-zero in the interval < ct < (1 -f ^7)^^, 
where it equals 



pia) = (NM^il - a)-'{a-' 1 h)'^'- (15) 
This leads to the effective number of degrees of freedom 



-4[(1 + 4/7)1/4 + (l-f 4/7) 



-1/41-2 



plotted in Fig. with a mean photocount of 



= 1^^/7(^1 + 477-1). 



(16) 



(17) 



For strong absorption, 7 ^ 1, we recover the black-body 
result i^cff = V, as expected. For weak absorption, 7 <C 1, 
we find — 2v^. In the weak-absorption regime we 
can compute the entire distribution P{n) analytically. 
The result 

P(n) (X (n"/n!)(l + /)-"/'if„-i/2 (^^off VlT/) , (18) 

with n = ^f\/j ^ ^ Bessel function, is Glauber's 
distribution [p_5| with a reduced number of degrees of 
freedom. 

Our second example is an optical cavity connected to 
a photodetector via an A^-mode waveguide. The cav- 
ity modes near frequency oj are broadened over a fre- 
quency range N lS.uj much greater than their spacing Aw 
ii N ^ 1. The cavity should have an irregular shape, or 
it should contain random scatterers — to ensure chaotic 
scattering of the radiation. For this system we define 
the normalized absorption rate as 7 = Tdweii/Ta, where 
Tdwcii = 2tt/NAuj ~ 1/wc is the mean dwell time of a 
photon in the cavity without absorption. The scattering- 
strength density for A^ ^ 1 follows from the general for- 
mulas of Ref. |2^. The result has a simple form in the 
limit 7 <C 1 of weak absorption. 



p{a) = {N/27r){l a)-\a - a) 



\l/2 



(19) 



for cr_ < (T < (T-)_ with cr± = 1 — 87 ± 2'^\/2. In the 
opposite limit 7 3> 1 of strong absorption, p{ij) is given 
by the same Eq. (^5|) as for the disordered slab. We find 
the effective number of degrees of freedom 

z/eff/i^= (1 + 7)^(7' + 27 + 2)-\ (20) 

plotted also in Fig. |^, with a mean photocount of 

n-z./7(l + 7)"'- (21) 

Again, VcS = for 7 3> 1. For 7 1 we now find 
VcS = \v. It is remarkable that the ratio t'oft/'^ for the 
chaotic cavity remains finite no matter how weak the ab- 
sorption, while this ratio goes to zero when 7 — > in the 
case of the disordered slab. 

These two examples concern thermal emission from ab- 
sorbing systems. As we discussed, our general formulas 
can also be applied to amplified spontaneous emission, 
by evaluating the Bose-Einstein function (|2|) at a neg- 
ative temperature. Complete population inversion cor- 
responds to / = — 1. A duality relation [E5l between 
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absorbing and amplifying systems greatly simplifies the 
calculation. The dielectric constants e' ± ie" of dual sys- 
tems are each others complex conjugates, so dual sys- 
tems have the same value of and 7. Their scattering 
matrices are related by Sj^ — S_ , hence the scattering 
strengths ci, (T2, . . . (Tat of an amplifying system are the 
reciprocal of those of the dual absorbing system. 

We need to stay below the laser threshold, in order 
to be in the regime of linear amplification. The semi- 
infinite medium is above the laser threshold no matter 
how weak the amplification p2| , but the cavity is below 
threshold as long as 7 < 1. We find that n and Vcs/^ 
are given by Eqs. ( |20| ) and ( |2l] ) upon substitution of 7 by 
—7. In Fig. H we compare Vfis/v for amplifying and ab- 
sorbing cavities. In the limit 7 — > the two results coin- 
cide, but the 7-dependence is strikingly different: While 
Vcs/i^ increases with 7 in the case of absorption, it de- 
creases in the case of amplification — vanishing at the 
laser threshold. Of course, close to the laser threshold 
[when 7^1 — (f^cTdwcu)^^^^] the approximation of a 
linear amplifier breaks down. 

In summary, we have derived a relation between the 
photocount distribution P{n), in the long-time limit, 
and the eigenvalues CTi , (72 , . . . (Jn of the scattering- matrix 
product S ■ S'^ . The super-Poissonian noise Varn = 
nil -\- n/vcs) is that of a black body with a reduced 
number v^s of degrees of freedom. We have computed 
VcS for several types of random media, in the large-A'^ 
regime, using results from random-matrix theory. In a 
weakly absorbing or amplifying chaotic cavity, the ra- 
tio Vcslv is a universal factor of 1/2 — independent of 
microscopic parameters. In a disordered slab, Vcsl^ van- 
ishes oc l/y^ for small absorption rates We have 
found that Vcs/v vanishes also on approaching the laser 
threshold in an amplifying chaotic cavity. 

The reduction of feff amounts to an excess noise of am- 
plified spontaneous emission. Its origin is the presence 
of a large number N of overlapping cavity modes, and 
a broad distribution p{a) of the corresponding scatter- 
ing strengths. Overlap of cavity modes is avoided in the 
usual laser geometry, but it is generic in a random laser. 
This fundamental difference was pointed out thirty years 
ago by Letokhov in the paper that pioneered the 
notion of a "stochastic resonator" . Letokhov concludes 
his paper by surmising that the statistical properties of 
spontaneous emission would be distinctly different from 
the usual case. The reduction of the number of degrees 
of freedom predicted here forms an experimentally acces- 
sible signature of this difference. 

I have benefitted from discussions with P. W. Brouwer, 
M. P. van Exter, and J. P. Woerdman. This work was 
supported by the Dutch Science Foundation NWO/FOM. 
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